Cryopreservation effect on DNA methylation profile in rainbow trout spermatozoa

Spermatozoa are the cells that are most commonly used for cryopreservation of valuable genetic resources in aquaculture. It is known that fish spermatozoa transmit to the embryo not only their genetic but also their epigenetic profile, especially DNA methylation. Therefore, any alteration of the DNA methylation profile in spermatozoa induces the risk of transmitting epigenetic alterations to the offspring. The aim of this study was to assess the effect of cryopreservation on DNA methylation in rainbow trout spermatozoa. To trigger variable cellular response after freezing–thawing, spermatozoa from mature males were cryopreserved with dimethyl sulfoxide, methanol or glycerol as cryoprotectant. We observed that dimethyl sulfoxide was the best to preserve thawed spermatozoa functions. Methanol only slightly preserved all the cellular parameters, while glycerol failed to protect motility and fertilization ability. The consequences on DNA methylation were assessed using Reduced Representation Bisulfite Sequencing (RRBS). Sperm cryopreservation did not thoroughly impact DNA methylation, although 335–564 differentially methylated cytosines were characterized depending on the cryoprotectant. Very few of them were shared between cryoprotectants, and no correlation with the extent of cellular damage was found. Our study showed that DNA methylation was only slightly altered after sperm cryopreservation, and this may render further analysis of the risk for the progeny very challenging.

and glycerol 11 .In all, cryopreservation might induce a risk of DNA methylation alteration in spermatozoa, but the extent of the damages at the scale of the genome and its functional regions has not been addressed yet.
Rainbow trout (Oncorhynchus mykiss) is important in aquaculture and it has been an ever-present species in selective breeding programs in the past decades 12,13 .Cryobanking of the improved commercial strains using cryopreserved sperm is operational in France 14 , but to date, the impact of sperm cryopreservation on DNA methylation has never been studied in rainbow trout.Therefore, testing the risk of DNA methylation changes at the spermatozoa level is the first step to be carried on before the consequences on the progeny can be considered.The aim of our study was then to assess whether cryopreservation of rainbow trout spermatozoa can alter their DNA methylation profile and whether alterations of the cellular quality of thawed sperm is related to possible DNA methylation alterations.To this end, we used different cryoprotectant molecules in order to induce some variability in cellular cryoprotection, namely DMSO, MeOH and glycerol.DMSO and MeOH are the most commonly used molecules for sperm cryopreservation in salmonids.However, DMSO is a highly chemically reactive molecule comparing to MeOH.Glycerol does not protect spermatozoa functionality as do DMSO and MeOH, but it still provides some structural protection.Thawed spermatozoa quality was evaluated from their membrane quality, mitochondria quality, motility and fertilization ability.DNA methylation was assessed at the genome-wide level using reduced representation bisulfite sequencing (RRBS), a method that allows a partial genome representation, thereby allowing the assessment of many biological replicates for each cryoprotectant.

Cellular quality of the thawed spermatozoa according to the cryoprotectant
Plasma membrane quality of the thawed spermatozoa was deduced from spermatozoa ability to exclude the membrane impermeant propidium iodide.As shown Fig. 1a, DMSO provided a good protection of the plasma membrane, since all samples had more than 75% of the cells with intact plasma membrane.Glycerol was also able to preserve the plasma membrane quality with an efficiency similar to that of DMSO.On the contrary, MeOH yielded the lowest percentages of cells with intact plasma membranes, and the highest variability in plasma membrane quality among the spermatozoa of the different males.Mitochondria quality of the spermatozoa was estimated from their JC-1 staining intensity (Fig. 1b).Again, DMSO and glycerol showed the best ability to preserve the mitochondria quality, with more than 50% of the cells having mitochondria with a high activity, while the spermatozoa cryopreserved with MeOH had a very low mitochondria quality.
Regarding motility, despite the overall low sperm motility after thawing (Fig. 1c), the highest values were obtained with DMSO.The spermatozoa cryopreserved in MeOH retained very little motility while glycerol completely fail to protect sperm motility ability.The same trend was observed with the fertilization ability of the cryopreserved spermatozoa (Fig. 1d).Our results showed that the best fertilization ability was obtained with DMSO.MeOH yielded lower performances, with no fertilization rates above 50%.After cryopreservation in glycerol, the spermatozoa completely lost their fertilization ability, which was to be expected when considering the total loss of motility of these samples.These fertilization results were confirmed in another fertilization test in which other straws from the same sperm samples were used, with another batch of oocytes (results not shown).
In all, cryopreservation with different cryoprotectants altered the functionality of the thawed spermatozoa, and the different cryoprotectants yielded different cellular response to cryopreservation.We observed that DMSO preserved the membrane quality and the mitochondria quality, and protected to some extent the motility and the fertilization ability.MeOH was less efficient at protecting whatever cellular quality level, while interestingly, glycerol was able to preserve the membrane quality and mitochondria quality, but was completely inefficient regarding the motility or the fertilization ability.This wide set of alterations depending on the samples provided us with a very variable biological material to explore further the risks of changes in DNA methylation after cryopreservation.

Validation of the RRBS sequencing data for DNA methylation analysis
Because RRBS yields a reduced representation of the analyzed genomes, we assessed whether our sequencing data still had enough cytosines in CpG sites to provide an acceptable representation of those in the whole rainbow trout genome.As a reference, we used a total number of cytosines (in CpG sites) of 70.6 × 10 6 , although the reference trout genome contains 35.3 × 10 6 CpG dinucleotides (Omyk_1.0genome version).Indeed, cytosines of both forward and reverse strands can be sequenced (named plus and minus strands depending on their sense on the reference genome), and this can result in the sequencing of the same CpG site twice.We did not remove this redundancy (+ and − CpGs) from our data, because we do not know whether a CpG site can be altered symmetrically by cryopreservation or not.As a consequence, the whole information was kept, and each sequenced cytosine (in a CpG site) in our conditions is one out of the 70.6 × 10 6 genomic cytosines (twice the number of genomic CpG sites).We found an average of 2.9 × 10 6 cytosines in CpG sites per sample (fresh and cryopreserved sperm), that represent 4% of the 70.6 × 10 6 genomic cytosines in CpG sites (Table1).This representation was slightly higher in MeOH, likely because these samples also had more reads than those in the other conditions, leading to more successful extraction of the cytosines.It is noteworthy that when combining all the cytosines of all the males from one condition, the representation of the genomic cytosines in CpG sites was doubled (9%).Hence, the reduced representation of the genomic cytosines in CpG sites by the data at the sample level was compensated by the large number of biological replicates in this study.
To investigate whether these 9% of genomic cytosines in CpG sites were homogeneously distributed over the entire reference genome, we analyzed their individual location.We observed a homogeneous distribution of the cytosines all along the chromosomes of the rainbow trout genome, with never more than 150 kb between two consecutive cytosines (Fig. 2).The same homogeneous distribution of the sequenced cytosines in CpGs sites was found within each sample of each cryopreservation condition.Therefore, this distribution indicates that www.nature.com/scientificreports/ We thus studied the individual methylation ratio of all the cytosines with a sequencing depth N ≥ 10.For all treatments, more than 80% of the analyzed cytosines in CpG sites had a methylation ratio of 1 (Fig. 3b, class 0.95-1).This confirms the highly methylated status of the cytosines in all samples, as shown above in the global results.The second most represented category encompassed the cytosines with a methylation ratio of 0. Only a very small proportion of cytosines had intermediate methylation values, meaning that there is very little variability in the methylation status of the cytosines between spermatozoa within a sample (most are in the binary status of either methylated or unmethylated).When we compared the fresh to the cryopreserved samples, we did not observe any shift between the 3 categories (Fig. 3b, percentages) that would have indicated that some cytosines would have acquired a more variable status among spermatozoa.Besides, no significant difference was found when comparing the 12 fresh samples to any of the cryopreserved ones (p > 0.05, Kruskal-Wallis test).Thus, if sperm cryopreservation has induced some subtle alterations of the cytosine methylation ratio, it was not detected from this representation.

Characterization of differentially methylated cytosines (DMCs) after sperm cryopreservation
The differential methylation of cytosines in CpG sites between fresh and cryopreserved samples was characterized using a paired-driven analysis, in order to focus on the cytosines within each male that were read in both the fresh and respectively the DMSO, MeOH and glycerol conditions.This prevented any statistical bias where neighboring cytosines are used in the comparison to compensate for the lack of information at a given site.These cytosines will be referred to as paired cytosines (i.e. a cytosine at a genomic position that is present in both the fresh sperm and the corresponding cryopreserved one).In this analysis, all cytosines with a sequencing depth N ≥ 1 were included.We first assessed how many paired cytosines would be shared between the different males.Surprisingly, the highest proportion of these cytosines (> 15%) were shared by 8-10 males (Fig. 3c) whereas fewer Table 1.Summary of the sequencing data obtained in the RRBS analysis.Total raw reads number includes reads 1 and reads 2 (paired end sequencing).Reads with adapter is the number of reads whose insert is shorter than 100 bp.The % mapping efficiency = (number of reads that mapped only once on the trout genome/ trimmed reads number) × 100.Sequencing depth is the number of times a cytosine has been sequenced, mapped to the reference genome and extracted as a cytosine in CpG sites.The % total genomic cytosines in CpG sites refers to (number of identified cytosines in CpG sites/number of cytosines in CpG sites in the reference genome-i.e.70.6 × 10 6 -) × 100.Values are expressed as mean ± SD of the 12 males in each condition (Fresh, Cryopreserved in DMSO, MeOH, or Glycerol).www.nature.com/scientificreports/cytosines (< 6%) were present in only 1-5 males.Furthermore, the number of paired cytosines shared from 6 to 12 males represented 90% of the total number of paired cytosines.That paired cytosines were found as being more abundant in more males would indicate that some fragments were more prone to be sequenced than others in the majority of the males, although the origin of this phenomenon regarding a possible sequencing bias is difficult to understand.In all, since we wanted to select the highest proportion of paired cytosines, we selected the paired cytosines present in at least 6 males (n ≥ 6).This threshold also allowed us to strengthen the biological meaning of the DMCs that we may obtain.When comparing DMSO, MeOH and glycerol to fresh sperm, we found 335, 564 and 357 DMCs, respectively (Fig. 3d).No DMCs were common to the three cryopreservation conditions, and only few DMCs were common to 2 cryoprotectants.None of these DMCs had a straightforward status as we observed that depending on the males, a given DMC could be hyper-or hypomethylated after cryopreservation.In order to sort out some trends in the methylation changes, the DMCs were classified as hypermethylated, or hypomethylated, when they had the same trend in > 50% of the males.As shown in Fig. 3e, cryopreservation with DMSO triggered the same number of hypomethylated and hypermethylated DMCs (35%) whereas in MeOH, more DMCs were hypomethylated (45%) than hypermethylated (33%).A reverse tendency was observed after cryopreservation with glycerol where more DMCs were hypermethylated (42%) than hypomethylated (32%).Some DMCs remained unclassified because they did not show a clear trend among the corresponding males.This emphasizes the fact that cryopreservation did not trigger a homogeneous effect in all the males, and that sperm cryopreservation impact on DNA methylation was different depending on the cryoprotectant.

Characterization of regions potentially sensitive to cryopreservation
Despite the low DMCs numbers obtained in our cryopreservation conditions, we investigated whether there were any DMRs (differentially methylated regions) between fresh and each cryoprotectant.No DMRs were found after cryopreservation with DMSO and MeOH, whereas only 2 DMRs were observed after cryopreservation with glycerol.Some potentially sensitive regions could have escaped detection by DSS because of the stringency of the chosen settings, where a DMR search requires that at least five consecutive CpGs are present in a sliding 50 bases window.In an empirical approach focused on the DMCs that we have identified, we sought for regions with at least two DMCs separated by no more than 100 bp (our maximum read size).These regions were referred to as potentially sensitive regions.We found 60 of these potentially sensitive regions, present in at least one comparison between cryopreserved sperm and fresh controls.Among them, 17 were sensitive to DMSO, i.e. they contained at least 2 DMCs in the DMSO-fresh comparison, 33 were sensitive to MeOH, and 16 were sensitive to glycerol.One region was affected by all 3 cryoprotectants, and 3 regions were sensitive to both DMSO and MeOH, but not to glycerol.
These regions found to be sensitive to the cryoprotectants were distributed in all genomic regions and these included features such as introns, exons, promoter 5kb (Table 2) and intergenic regions (33 regions, Supplementary Table S1).Some annotated genes were identified in these regions, although the only region common to the 3 cryoprotectants belonged to an intergenic region (8: 79,446,507-591).It is noteworthy that among the regions related to annotated genes (promoters, introns and exons), only 4 of them were affected by DMSO while 12 were affected by MeOH and 12 by glycerol.We found that the distal promoters of ebf2 and p4ha1a genes were affected by cryopreservation with DMSO.Among the promoter regions affected by MeOH, one regulates u2 gene.The other affected distal promoters were localized in the mitochondrial chromosome.These promoters are related to cox2, cox3, nd3, nd4, nd4l, atp6 and atp8 genes.The pcdhb and slc26a2 genes had their first intron affected as well.We also found that among the distal promoters affected by glycerol, one was associated to znf385d gene and another one was associated to osbpl10 gene.Other genes such as ptprub and ndufs8a were also affected at their introns by glycerol.

Relationship between DNA methylation changes and cellular alterations
To explore the possible link between DNA methylation changes after cryopreservation and the corresponding cellular alterations of the same samples, we assessed whether there was a correlation between the number of DMCs and the fertilization rate.Indeed, fertilization ability is a cellular quality parameter that combines several functions of the spermatozoa.This includes its motility capacity, its ability to reach the egg micropyle, its plasma membrane ability to fuse with the egg membrane, its nucleus ability to decondense and to fuse with the female pronucleus, and finally the paternal DNA ability to sustain development up to eyed stage, far later than embryonic genome activation.No significant correlation (p > 0.05) was found between the 2 parameters.Because glycerol did not yield any fertilization, thus preventing the testing of some correlation, two other cellular parameters that are plasma membrane quality and mitochondrial activity were also tested.No correlation with DMCs number was found either, whatever the cryoprotectant.Surprisingly, a significant negative correlation was found between DMC number and motility after cryopreservation with DMSO (p = 0.04, not shown).

Discussion
Cryopreservation of rainbow trout sperm with 3 different cryoprotectants led to varying pattern of cellular alterations of thawed spermatozoa in our study.Although cryoprotectant comparison for rainbow trout sperm cryopreservation was not the scope of this study, the hierarchy between cryoprotectants that we observed in our study deserves some comments, especially regarding the poor performances of MeOH and glycerol.Indeed, fertilization rates (at eyed stage or at hatching) above 80% have been reported after cryopreservation with MeOH [15][16][17] , yielding performances comparable to those with DMSO [18][19][20] .In their thorough comparison, Babiak et al. 21were proposing that dimethyl acetamide was yielding the best fertilization rates, followed by glycerol, then methanol, but they were also showing that the best cryoprotectant could lose its rank depending on the additives.www.nature.com/scientificreports/It is also known from a model study that the extender composition can reverse the cryoprotectant efficiency 22 .
In a previous work, we had observed that MeOH was yielding higher fertilization rates than DMSO 23 , but the additives were different (egg yolk instead of liposomes), and we have never been able to repeat this results with the conditions of the present work.Why MeOH yielded such low sperm motility and fertilization performances after thawing in our study might indeed be due to the composition of our extender (containing liposomes and sucrose).Regarding the plasma membrane integrity parameter, MeOH 24 and DMSO 25,26 are known to alter plasma membrane fluidity and permeability.This is why the spermatozoa were washed before the tests, but it cannot be excluded that propidium iodide exclusion was lessened in membranes previously exposed to MeOH, and this would misleadingly increase the percentage of cells with altered membranes.However, it is difficult to ascertains why the sperm plasma membranes would be more affected by the remnants of MeOH than by those of DMSO, as both are short acyl chain molecules.Despite our inability to clearly understand the roots of the differences between the tested cryoprotectants, these differences validated our choice of the cryoprotectants, as they induced different cellular responses to cryopreservation.In addition, the variable response towards cryopreservation among fish males revealed an intrinsic biological variability.Such overall variability provided us with a set of samples that was suitable to study the effect of cryopreservation on the molecular parameter that is DNA methylation.The RRBS method used in this study yielded a limited number of cytosines in CpG sites, and this reduced set of data made possible the use Table 2. List of the potentially sensitive regions according to the cryoprotectant.The sequenced regions which displayed at least 2 DMCs (differentially methylated cytosines) between fresh and cryopreserved samples within 100 b (read size) in regions others than intergenic were included in this list.Gly: glycerol.Regions: chromosome number and region coordinates on the reference genome.Gene: name from the ensembl database (version 104) according to the reference genome (Omyk_1.0genome version).Promoter 5 kb: region 5 kb upstream of the transcription starting site; Exon 1: first exon; Downstream: region 1 kb downstream of the transcription end site; NA: not available.*the 2 DMCs belonged to the exon1 of cox1, downstream of nd2, promoter 5 kb of cox2, promoter 5kb of atp8, promoter 5kb of atp6, promoter 5kb of cox3, promoter 5kb of nd3, promoter 5kb of nd4l, promoter 5kb of nd4.www.nature.com/scientificreports/ of many biological replicates to explore the impact of sperm cryopreservation on DNA methylation.Our results demonstrated that cryopreservation had no straightforward impact on the global DNA methylation, although it yielded some DMCs between cryopreserved and fresh sperm.Presumed regions that would be more sensitive to cryopreservation were proposed, but few were shared between cryoprotectants.
Our study showed that the extent of the cellular cryoprotection provided by the different cryoprotectants had no related consequence on DNA methylation alteration.Indeed, using either an optimal (DMSO), non-optimal (MeOH) or the worst cryoprotectant (glycerol) did not affect the global DNA methylation in rainbow trout.This is at odds with some other species such as goldfish 10 and European eel 9 , where the optimal cryoprotectant which is MeOH did not alter global DNA methylation assessed by LUMA, while DMSO induced a hypomethylation.Therefore, our study revealed the stability of global DNA methylation towards the cryopreservation whatever the cryoprotectant used in rainbow trout.However, since we used RRBS technique, this allowed us to go further than global methylation and to assess DNA methylation at the nucleotide level.
We observed a limited number of DMCs after sperm cryopreservation, although a high number of paired cytosines was analyzed in each comparison.Besides, we were not able to find any DMR except 2 in glycerol.This further emphasizes how slight is the impact of cryopreservation on sperm DNA methylation in rainbow trout.This failure to detect more DMCs and DMRs is not due to the reduced representation of RRBS used here.Indeed, in a study dealing with rainbow trout physiological adaptation to its environment, 108 DMRs were found in spermatozoa between hatchery and natural origin rainbow trout using the same RRBS technique 27 .We surmise that the impact of hatchery on DNA methylation was important enough to induce many DMCs leading to many DMRs.In our study, since the impact of sperm cryopreservation is tenuous, it yielded a restricted number of DMCs and as a consequence the DMCs were not dense enough to form DMRs.However, despite the limited number of DMCs found in our study, these DMCs are reproducible among the males since they were found in at least 6 males.This reveals the importance of having many biological replicates when assessing the impact of cryopreservation in order to validate the slight impact that can be observed.
Our analysis revealed that cryopreservation induced both hypomethylation and hypermethylation depending on the DMCs.As a consequence, no clear trend could be proposed as to the effect of cryopreservation.In addition, the status of the same DMC differs among the males, and the same male has DMCs with different status.Because of the quiescent state of sperm chromatin, the methylation changes were unlikely mediated by the enzymes responsible for methylation (the DNMTs) and demethylation (the TET pathway).Besides, the process of freezing-thawing of the spermatozoa is too fast for any enzyme being able to catalyze the methylation or demethylation reaction.The hypomethylations and hypermethylations observed in our study were more likely resulting from chemical impact on DNA methylation.Until now, it is still unclear how the cryoprotectants are able to modify the DNA methylation.DMSO is a chemically reactive molecule bearing 2 methyl groups, and it has been shown in an in vitro study that DNA can acquire methyl groups when incubated with DMSO and hydroxyl radicals 28 .In the case of cryopreservation, the releasing of reactive oxygen species by damaged sperm mitochondria may trigger the formation of these hydroxyl radicals and induce this DMSO-associated methylation event.
Regarding the demethylation event, it has been shown that 5-mC can be oxidized into 5-hmC in the presence of reactive oxygen species 29,30 .The further transformation of 5-hmC into cytosine was shown to be made possible at high pH 31 .It is known that pH can increase upon freezing of buffered solutions 32 , and it cannot be excluded that a similar change is occurring locally within the cryopreserved spermatozoa.Therefore, demethylation may take place by chemically oxidizing 5-mC, and 5-hmC transformation into a cytosine may be linked to local pH changes during cryopreservation.The puzzling part in all these chemically-induced methylation changes is that both methylation and demethylation took place in the same sperm sample, and we have no clue as how such microreactors are made possible, either within the same nucleus or between nuclei.
Despite the lack of DMRs in our study, we could identify potential sensitive regions in which several cytosines showed methylation changes in one or two cryoprotectants.Some of these regions are related to gene promoter regions that are instrumental to gene proper expression pattern.We focused our discussion on the genes affected by DMSO and MeOH, as glycerol cannot produce any offspring.Interestingly, the genes of promoter regions affected by cryopreservation with DMSO (ebf2, p4ha1a) are involved in specific processes of embryogenesis.In zebrafish, ebf2 (early B cell factor 2) is a transcription factor involved in primary neurogenesis in the embryo 33 , where it stabilizes the commitment of the neuronal progenitors.The p4ha1a gene (prolyl 4-hydroxylase, alpha polypeptide I orthologue a) is an ortholog of p4ha1 whose protein is involved in the regulation of early growth in bone tissues in bighead carp 34 .Regarding the genes whose introns or promoters were affected after cryopreservation with MeOH, some had developmental functions as well (pcdhb, slc26a2, u2).While no function of pcdhb and slc26a2 was described yet in fish, it has been shown in mice that pcdhb is involved in the regulation of neural development during the axonal assembly 35 and the axonal convergence into glomeruli in the olfactory bulb 36 .Expression inhibition of slc26a2 (diastrophic dysplasia sulphate transporter) in mice is altering bone growth 37 , and aberrant methylation of slc26a2 in human inhibits this gene expression and reduces its protein expression 38 .Concerning u2, it has been shown in zebrafish that its protein belongs to the pre-mRNA Retention and Splicing (RES) complex that is fundamental during zebrafish embryogenesis 39 .The other genes (cox2, cox3, nd3, nd4, nd4l, atp6 and atp8) are involved in the oxidative phosphorylation (OXPHOS) reaction in fish as they belong to different mitochondria respiratory complexes [40][41][42] .In light of the limited amount of data that was collected for each region, these affected genes represent preliminary clues about a slight impact of cryopreservation on DNA methylation, and provide a basis for further studies on important candidate genes.Recently, one study claimed that differentially methylated genes would be found after cryopreservation of Black Rockfish sperm 43 , but the authors studied CpC dinucleotides instead of CpGs.These dubious methylation results on sites which are not supposed to be methylated in vertebrates unfortunately prevent any comparison with our present data.Indeed, our study was focused on the CpG sites, the only sites whose changes in methylation status are known to have functional consequences in vertebrates.www.nature.com/scientificreports/Our data revealed the absence of a link between alteration of DNA methylation and fertilization ability.The males that had the worst quality for a cellular quality parameter did not necessarily have the highest number of DMCs.This indicates that the damages which altered the spermatozoa functions to a great extent, such as cell dehydration and ice crystal formation 44 , had little effect at the chromatin level.As a consequence, none of the cellular parameters used in this study can be used as markers of the risk of DNA methylation alterations.This lack of correlation can be explained by the general stability of the DNA methylation, validated by the low number of DMCs and the overall maintenance of the global DNA methylation level in all samples.
At this point, it is difficult to estimate whether the slight alterations observed in our study bear any consequence for the progeny.Indeed, our sequencing data reveal the DNA methylation pattern of the whole spermatozoa population, whereas fertilization selects a minor proportion in this population.Among the thawed spermatozoa, there may be some whose DNA methylation was not affected by cryopreservation.If one of these spermatozoa is the one which fertilizes an oocyte, there will be no transmission of any DNA methylation alteration.Nevertheless, if a spermatozoon with an altered DNA methylation pattern did fertilize an oocyte, there may be a risk of transmitting this alteration to the progeny, although the extent of such transmission cannot be predicted.It is known that DNA strand breaks can be repaired by the maternal machinery in rainbow trout, and that in some cases, inhibiting the repair mechanism will increase the post-hatching mortality rate 45 .However, nothing is known yet about DNA methylation repair after fertilization, and what would be the template to direct this repair.It was shown in zebrafish that specific histone variants in the spermatozoa (H2AFV) transmit the hypomethylated status of some DNA loci to the early embryo 3 .These variants stabilize the DNA hypomethylated status of these loci by deterring DNA methylation.However, we do not know whether these variants can protect their DNA loci from the spurious methylation induced by cryopreservation.If not, these variants may still be able to redirect a proper DNA demethylation in the embryo, as they did in maternal haploids 3 .The existence of such a speculative DNA methylation repair mechanism in the embryo would be strengthened if demonstrations were made that spermatozoa with experimentally altered DNA methylation do not sire abnormal or unviable embryos.However, to the best of our knowledge, we could not find such evidences, apart from the study on xenopus embryos sired with spermatids, whose DNA methylation profile differs from that of the spermatozoa, and where abnormal and reduced developments were observed 46 .
Regarding the slight modifications obtained after cryopreservation in our study, if the offspring is to inherit the altered DNA methylation pattern of the spermatozoa due to the lack of repair mechanism, and that the affected genes are essential for embryo development, this should induce some increased embryonic mortalities.However, fertilization with cryopreserved sperm in rainbow trout was reported to be harmless regarding embryo development, once fertilization is triggered 47,48 .It cannot be excluded however that non-lethal long-term effect could be transmitted.For example, rainbow trout embryos sired with cryopreserved spermatozoa had altered expression of several marker genes 49 , reduced growth two weeks after hatching 50 , and their stress resistance was modified at later stages 51 .In all, the effect of sperm DNA methylation changes on the progeny is utterly difficult to assess.Indeed, when the effect is not immediately made visible by increased embryonic mortalities, or by obviously reduced growth rates after hatching, the biological targets can be multiparametric, or adopt a wide range of symptoms in a wide range of tissues over the whole life of the individual.This is the case for example with artificial reproductive technologies (ART) in human (reviewed in 52 ), including when in some cases, ARTinduced epigenetic defects resolves at later age 53 .In another context, changes in zebrafish sperm DNA methylation (induced by F0 male exposure to mercury) were associated with far off symptoms that were hyperactivity and retinal dysfunction of their adult progeny (F2) 54 .Farmers may be more interested in checking stress and disease coping abilities, feed conversion and growth ability etc., but such criteria, already assessed in selective breeding schemes, are often very difficult to separate from the environmental specificities the individual were reared in 55 .Besides, it is now accepted that spermatozoa is transmitting a large variety of information to the progeny, including clues on the life of the father (reviewed in 56 ), thereby adding another difficulty to assess the extent of the transmission that may not be solely related to DNA methylation alteration.Additionally, it was shown in Atlantic salmon that small changes in the fertilization procedure, for example spermatozoa number, induced some changes in the offspring characteristics 57 , likely by changing the fertilizing sub-population of spermatozoa.Another unexpected difficulty to sort out epigenetic effect was shown in cattle, where sperm DNA was experimentally damaged (oxidative stress) 58 .The authors showed that, although the damage had no link with DNA methylation changes in the spermatozoon, it reduced the paternal chromatin ability to be demethylated during early development of the embryo.
In all, as long as epidemiological studies are not set up, it is only when altered genes are identified in the gametes that a candidate analysis can be performed in the progeny.This enforces the strategy chosen in the present work, that was to test in the spermatozoa first whether some important genomic sites would be repeatedly and consistently affected by cryopreservation.We observed only marginal modifications, that were variable in their expression between individuals, and that were not related to the intensity of the cellular damages induced by cryopreservation.At this stage of our knowledge, this would indicate that cryopreservation of rainbow trout spermatozoa is not triggering DNA methylation damages that would bear a high risk for the progeny.

Ethics statement
Rainbow trout males from the synthetic strain were kept in the experimental fish facilities ISC LPGP of INRAE (Agreement number D-35-238-6) with full approval for experimental fish rearing in strict line with French and European Union Directive 2010/63/EU for animal experimentation.The Institutional Animal Care and Use Ethical Committee in Rennes LPGP (Fish Physiology and Genomics Department) specifically approved this study (n° T-2020-37-CL).All fish were handled for gamete collection in strict accordance with the guidelines Vol:.( 1234567890

Sperm collection and cryopreservation
Rainbow trout from the SYNTHETIC INRAE strain were sired in January 2019 at the experimental fish farm UE PEIMA of INRAE.After development in natural spring water in the dark up to hatching, they were transferred into indoor tanks supplied with the same water, under natural photoperiod.At 4 months post fertilization, they were transferred into outdoor tanks supplied with dam water.They were fed ad libitum with commercial pellets (Le Gouessant B-Repro 32 ASTX-semi F7).After 2 years (January 2021) the males (average weight 1.2 kg) were transferred to the ISC LPGP fish facility.After 2 weeks in recirculated water tanks (fiberglass, 12 °C, winter photoperiod, same feeding as in the fish farm), 12 fish were anaesthetized, and their sperm was collected by gentle striping and stored individually on ice.All sperm samples were treated individually (12 biological replicates in every analyses).They all had high motility and high membrane quality percentages (mean of the 12 males: 99.92 ± 0.2 and 98.05 ± 0.84 respectively).Fresh sperm concentration was adjusted to 6 × 10 was snap-frozen in liquid nitrogen prior to storage at − 80 °C.This fraction was referred to as fresh sample, as it was not submitted to the whole cryopreservation process (no exposure to cryoprotectant, no slow freezing and thawing and no associated cellular constraints 44 ).The remaining sperm sample was split into 3 fractions to be cryopreserved with one of the three penetrating cryoprotectants dimethylsulfoxide (DMSO), methanol (MeOH) and glycerol (Fig. 4).Sperm cryopreservation was carried out with the FreezeFish kit (IMV Technologies ref 026913).One volume of sperm at 6 × 10 9 cells/ml was mixed with one volume of each cryopreservation solution (buffered saline with sucrose and liposomes and one of the three penetrating cryoprotectant: DMSO, MeOH or glycerol).The final cryoprotectant concentration with the spermatozoa was 7.5% (vol %) and the final concentration of liposomes (Freezelip) was 1.75% (vol %).The mix of sperm suspension with each cryopreservation solution was loaded into 500 µL French straws (36-48 straws per male and 144-192 per condition).Straws were cryopreserved in the programmable freezer Micro DigitCool according to a fast-freezing curve (bovine curve, IMV Technology, − 4°/min up to -12 °C, −40 °C/min up to -110 °C, and − 20 °C/min up to − 140 °C) and were then stored into liquid nitrogen until use.All males were treated at the same time in the same cryopreservation session.Thus, the time between sperm collection and sperm treatment (60 min), and between sperm mixing with each of the 3 cryoprotectant and cryopreservation (40 min) was the same for all males.www.nature.com/scientificreports/When needed, the straws were thawed in a water bath at 37 °C for 10 s.Thawed sperm was then washed by dilution with SFMM (1 vol + 20 vol) and centrifugation at 200 g for 10 min at 4 °C.The pellet was resuspended at 100 × 10 6 cells/mL in SFMM unless otherwise stated, and maintained on ice before use.

Membrane quality
Membrane integrity of the spermatozoa was assessed from their ability to exclude the fluorescent probe propidium iodide.Fresh and washed frozen-thawed sperm samples diluted at 1 × 10 6 cells/mL were incubated with propidium iodide (final concentration 12 µM) for 2 min and the percentage of unstained cells (cells with intact plasma membranes) was assessed by flow cytometry (MACSQuant Analyzer 10, Ref 130-096-343).Measurements were performed in duplicate on each sperm sample, using two different sperm suspensions in propidium iodide, within 1 h after thawing.

Mitochondria quality
Mitochondrial membrane potential of the sperm samples was assessed with the JC-1 fluorescent probe using the kit "MitoProbe™ JC-1 Assay Kit for Flow Cytometry (M34152, Invitrogen by Thermo Fisher Scientific).Diluted sperm in SFMM was incubated for 15 min at 20 °C with JC-1 at a final concentration of 0.5 µM.JC-1 fluorescence was then measured by flow cytometry (MACSQuant Analyzer 10).The population with the highest red and green fluorescence was selected as the one whose spermatozoa had high quality mitochondria, and its representation was given as a percentage of the whole sperm population.Each sperm sample was measured in triplicates within 1 h after thawing.

Motility
Thawed spermatozoa motility was assessed using a Computer-assisted Sperm Analysis (CASA) IVOS II using the software HTCasaII (IMV Technologies).The measurements were performed in triplicates for each sperm sample within 3 min post thawing on unwashed spermatozoa.For each replicate, 10 µL of thawed sperm (2.4 × 10 9 spz/mL) were activated with 600 µL of Actifish (IMV Technologies), and a fraction (2.8 µL) was loaded within 10 s into the chamber of a Leja slide (Leja®, REF.025107-025108).Motility was recorded from 10 to 15 s post activation (8 movies of 30 pictures at a speed of 60 pictures/min), and spermatozoa with an average path velocity (VAP) ≥ 30 µm/s were counted as motile.Results are expressed as a percentage of motile spermatozoa.

Fertilization ability
The fertilization test was conducted by adding 100 µL of thawed unwashed spermatozoa and 10 mL of Actifish (IMV Technologies) to approximately 150 oocytes (2 × 10 6 spermatozoa/oocyte).As a control for the quality of the oocytes, 25 µL of fresh sperm and 10 mL of Actifish (IMV Technologies) were added to approximately 150 oocytes every 8 fertilization tests.All sperm samples (12 males × 3 cryopreservation conditions and the fresh controls) were tested for fertilization on the same egg batch obtained from the pool of 30 spawns (spring strain females from UE PEIMA, 3 years old).The spawns were collected one day before the experiment and transported over night at 4 °C under oxygen to ISC LPGP.After fertilization, eggs were incubated in the dark at 10 °C in recirculated water at the fish facility ISC LPGP.Fertilization rates were assessed from the development rates at eyed stage (about 20 days post fertilization).To do so, the eggs of each fertilization test were transferred into a petri dish and a picture was taken.The total number of eggs was counted with the macro VisEgg 59 on Fiji-ImageJ while the number of fertilized oocytes (embryos with developed eyes) was manually counted with the Fiji-ImageJ counter.Fertilization percentage was calculated by the following formula: (number of fertilized eggs/total number of eggs × 100).The normalized fertilization percentage used in the results correspond to: fertilization rate of each sample/average fertilization rate of fresh sperm × 100.

Statistical analyses
Significant differences between cryoprotectants were assessed with the non-parametric Kruskall and Wallis test (multiparametric) followed by a two ways comparison with the Mann and Whitney test.Significance threshold are indicated in the figures, but the highest p value accepted is 0.05.Data representation was made using the non-parametric whisker box-plots, the whiskers indicating the minimum and maximum values, and the box showing, from bottom to top, the first quartile, the median, and the third quartile.

DNA extraction
Fresh spermatozoa previously snap-frozen (50 × 10 6 spermatozoza) were thawed by adding 800 µL of warm (42 °C) DNA releasing medium containing TNES buffer (125 mM NaCl, 10 mM EDTA, 17 mM SDS, 4 M urea, 10 mM Tris-HCl, pH 8) and proteinase k (P6556; Sigma-Aldrich) at 75 µg/mL (final concentration).Thawed washed spermatozoa (4 × 10 6 spermatozoa) were mixed with 800 µL of DNA releasing medium.All samples were digested overnight at 42 °C.After cooling, 800 µL of phenol:chloroform:isoamyl alcohol (25:24:1 vol) were added to the 800 µL of digested spermatozoa.After gentle shaking for 15 min at room temperature, the two phases were separated by centrifugation for 15 min at 8000g at 4 °C.A fraction of 400 µL of the upper aqueous phase was mixed with 100 µL NaCl 5 M in water and 1 mL cold 100% ethanol and incubated for 15 min at room temperature prior to centrifugation for 15 min at 12,000g at 4 °C.The DNA pellet was washed with 100 µL 75% ethanol and centrifuged for 15 min at 12,000g at 4 °C.The new DNA pellet was left to air dry.It was then resuspended in RNAse (A7973; Promega) 100 µg/mL final concentration in nuclease-free water and incubated for 1 h at 37 °C.DNA concentration was measured with the Qubit3 Fluorometer-invitrogen by Thermo Fisher Scientific using the Qubit™dsDNA BR Assay kit (REF Q32853, invitrogen by Thermo Fisher Scientific) according to the manufacturer's instructions.
For practical reasons, the same thawed sperm was used for membrane quality assessment and for DNA extraction.In parallel, the same thawed sperm but from different straws was used for motility and mitochondria quality assessment.

RRBS library preparation, sequencing and data processing
DNA digestion with MspI, library preparation and bisulfite conversion of the 48 samples (12 males × 4 conditions) using Diagenode Premium RRBS kit were performed by the Integragen company (Evry, France).Paired-end sequencing (100 b read length) was carried out by the company on Illumina NovaSeq™ 6000, with S2 patterned flow cells.Because of some difficulties with the first sequencing round (too many reads with adapters, indicating too short inserts), some samples have been resequenced after some changes in fragment size selection.These yielded an average of 66.5 × 10 6 paired end raw reads (133 × 10 6 total raw reads).The sequencing quality of the raw data was checked with FastQC program (version: FastQC_v0.11.7) (https:// www.bioin forma tics.babra ham.ac.uk/ proje cts/ fastqc/).Raw reads were trimmed using Trim Galore to remove adaptors sequences, the bases added during RRBS library preparation, and the bases with a phred score < 30 (https:// github.com/ Felix Krueg er/ TrimG alore version: TrimGalore-0.6.5).In order to exclude any risk of overestimation of DNA methylation due to a low bisulfite conversion (< 98%), we calculated the bisulfite conversion rate and found a mean value 99.6%.
Using bowtie2 aligner (version: bowtie 2-2.3.5.1) available in Bismark bio-informatic software 60 , trimmed reads were aligned to the reference genome (Omyk_1.0genome version).Methylation status of every cytosine in CpG sites was extracted using "bismark_methylation_extractor".These steps are available in the workflow which agrees to FAIR principles and is accessible online (https:// forge mia.inra.fr/ lpgp/ methy lome).All the samples had a median sequencing depth of the cytosines in CpG sites above 18 (Table 1) indicating that the data per sample was sufficient.To visualize methylation levels distribution on the reference genome, BigWig files were obtained with bedGraphToBigWig 61 and methylation values were visualized with Integrative Genome Viewer 62 .Because both 5-mC and 5-hmC are protected from bisulfite conversion into uracil 63,64 , the cytosines claimed to be methylated in this RRBS study are in fact either methylated or hydroxy methylated, and the information about the proportion of 5-hmC is not available.However, because 5-hmC are mainly found at an early step of the demethylation process (reviewed by Kohli and colleagues 65 ), we consider that their merging with 5-mC did not induce any bias in the interpretation of the data.
For the identification of differentially methylated cytosines (DMCs), the paired design of DSS general experimental design from DSS 66,67 was used.For the DMCs detection, all cytosines with a sequencing depth N ≥ 1 were included in the analysis.Indeed, almost all cytosine sites are either completely methylated or unmethylated (Fig. 3b), meaning that averaging the methylation status of several cytosines at the same site is not absolutely mandatory.In addition, DSS was also used to detect the differentially methylated regions (DMRs).The criteria used to determine the DMRs were to have a sliding sequence of at least 50 bases with ≥ 5 CpGs with a minimum of 75% DMCs (with p-value < 0.001).Ven diagram was created using Jvenn 68 , a web-based tool.The position of the DMCs within genomic features and the identity of the genes corresponding to regions potentially sensitive to cryopreservation was obtained by annotating DMCs with GenomeFeatures 69 and ensembl v104.1 annotation.

Figure 1 .
Figure 1.Cellular quality of the thawed spermatozoa according to the cryoprotectant.(a) Cells with an intact plasma membrane, analyzed from the percentage of cells which excluded the impermeant fluorochrome propidium iodide; (b) Cells with high mitochondrial activity, analyzed from the percentage of cells which had the highest fluorescence intensity after JC-1 labelling; (c) Motile cells calculated from the percentage of cells which had an average path velocity ≥ 30 µm/s; (d) Development at eyed stage after fertilization with cryopreserved spermatozoa.The percentages were normalized according to the development rate of the egg quality controls (fertilized with standard fresh sperm).In the boxplot representation (a-d), each dot corresponds to individual value for one sperm (n = 12 biological replicates in each condition).The boxplots show (from the bottom to the top): minimum, 25% percentile, median, 75% percentile and maximum.Significant differences between cryoprotectants (Kruskall and Wallis followed by Mann-Whitney test) are indicated by *p value < 0.05; **p value < 0.01; ***p value < 0.001, ****p value < 0.0001.

Figure 2 . 2 Figure 3 .
Figure 2. Example of the homogeneous distribution of the sequenced cytosines in CpG sites over the reference genome.The displayed region is representative of the pattern observed with all the sequenced data over the whole genome.The chart was obtained with IGV_2.8.13 and shows the distribution on the chromosome 1 from the position 58,234,274-59,476,270 (1.2 × 10 6 bases) of the rainbow trout reference genome "Omyk_1.0genome version".Each line corresponds to a different condition in which all the sequenced cytosines (in CpG sites) of all 12 males are represented.Each sequenced cytosine in a CpG site is represented by a bar.The bottom line represents the genes (dark blue) associated to this genomic fraction.

Figure 4 .
Figure 4. Experimental design.One sperm sample from each rainbow trout male (n = 12) was split into 4 fractions: one fraction was snap frozen in liquid nitrogen and used as fresh control.The 3 other fractions were cryopreserved, each one with a different cryoprotectant: dimethyl sulfoxide (DMSO), methanol (MeOH), and glycerol.

Condition Total raw reads × 10 6 Reads with adapter × 10 6 (% of total) % mapping efficiency Median sequencing depth of cytosines Number of identified cytosines × 10 6 % total genomic cytosines in CpG sites
www.nature.com/scientificreports/ of the Institutional Animal Care and Use Ethical Committee in Rennes LPGP (Fish Physiology and Genomics Department).CL is accreted by the French Veterinary Authority for fish experimentation (no.005239).The animal study is reported in accordance with ARRIVE guidelines (https:// arriv eguid elines.org) for animal research.